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We provide a theoretical study of the conductance response of systems based on graphene nanorib¬ 
bon to the potential of a scanning probe. The study is based on the Landauer approach for the 
tight-binding Hamiltonian with an implementation of the quantum transmitting boundary method 
and covers homogenous nanoribbons, their asymmetric narrowing and quantum point contacts of 
various profiles. The response maps at low Fermi energies resolve formation of n-p junctions induced 
by the probe potential and a presence of zigzag-armchair segments of the edges for inhomogeneous 
ribbons. For an asymmetric narrowing of the nanoribbons the scanning probe resolves formation 
of standing waves related to backscattering within the highest subband of the narrower part of the 
system. The QPCs containing a long constriction support formation of localized resonances which 
induce a system of conductance peaks that is reentrant in the Fermi energy, with the form of the 
probability density that can be resolved by the conductance mapping. For shorter constrictions 
the probe induces smooth conductance minima within the constrictions. In general, besides the 
low-energy transport gap, in the wider parts of the ribbon the variation of the conductance of the 
system is low compared to the narrower part. 


I. INTRODUCTION 

The electron flow in semiconductor nanostructures 
containing a two-dimensional electron gas (2DEG) buried 
shallow beneath the surface of the structure can be 
probed by a charged tip of an atomic force microscope. 
The technique, known as the scanning gate microscopy 1 
(SGM), gathers conductance maps as functions of the po¬ 
sition of the tip that is capacitively coupled to the elec¬ 
tron gas. For 2DEG confined in bulk semiconductors the 
maps are used to extract electron trajectories, branching 
of the electron flow via quantum point contacts (QPCs)®!, 
many-body phenomengP! Aharonov-Bohm effectS and 
scarred wave functions in quantum billiardS as well as 
the interference effects' 3 due to backscattering by the tip. 
The SGM technique has also been usecf3ED f or graphene 
with the electron gas remaining strictly on the surface of 
the structure. In particular, SGM technique was applied 
in studies of charge inhomogeneitieS^, charging local¬ 
ized states formed within a widening of the nanoribbonS 
and at potential constriction^®, universal conductance 
fluctuation^ 11 ! and weak localization effects 12 ! Very re¬ 
cently a SGM study of the graphene quantum rings was 
reported® with conductance fringes as due to resonant 
localized states. 

Quantum point contacts are basic elements of the 
quantum transport circuitry allowing for the current in¬ 
jection and detection, control of the number of con¬ 
ducting of modes, etc. The conductance quantization 
for QPCs defined within the bulk semiconductors is 
usually observed for electrostatic potentiald 14 ! control¬ 
ling the width of the constriction. For graphene sys¬ 
tems this approach is rather excluded, since the Klein 
tunneling 15 16 and perfectly conducting channel! 17 * 1 in 
graphene nanoribbons lead to a low effectiveness of elec¬ 
trostatic confinement and backscattering, allowing the 


Dirac electrons to pass across potential barriers 1 -®!. The 
QPCs are therefore for med by t ailoring constriction^®!® 
of the nanoribbond®^ 27 ^®®!^. The ineffectiveness of 
backscattering in graphene seems reflected by the results 
of the experime ntal SGM studies of the quantum point 
contacts (QPC]f®®®! which find flat conductance maps 
outside the QPC® in a distinct contrast to the results 
obtained for bulk semiconductor in which interference 
fringes due to the backscattering by the tip and the re¬ 
sulting formation of standing wave between the QPC and 
the probe are clearly observecP. The scanning gate mi¬ 
croscopy was used 5 ®! f or detection of spontaneous quan¬ 
tum dots formed along the disordered ribbon within the 
transport gap—. 

The purpose of the present paper is to determine the 
response of graphene quantum point contacts formed by 
constrictions of the nanoribbons to the scanning probe 
[Fig. 1(c)]. As a starting point of the study we consider 
reaction of homogenous graphene nanoribbons [Fig. 1(a)] 
to the perturbation by an external potential and asym¬ 
metric narrowing [Fig. 1 (b)] of the channels!®. The 
QPC [Fig. 1(c)] is formed by two asymmetric connections 
[Fig. 1(b)], whose scattering properties and response to 
the probe needs to be determined in order to separate the 
effects of the resonances formed within the QPC. We con¬ 
sider systems of perfect armchair or zigzag edges - and 
systems varying with atomic-step edges, as can be pro¬ 
duced with the Joule heating technique 35 . For pristine 
graphene ribbons®®!® the edge determines the charac¬ 
ter of the ribbon (semiconducting or semimetallic) and 
the ribbons with boundaries other than armchair have 
the transport properties similar to the zigzag edged®. 
For the purpose of the present study we develop an im¬ 
plementation of the quantum transmitting boundarjf 37 ! 
method for the solution of the coherent transport prob¬ 
lem set by scattering solution of the Schrodinger equation 
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for the tight-binding Hamiltonian. We study both the ef¬ 
fects of the local potential variation induced by the tip 
on edge-localized states and the effects of backscattering 
within the constrictions. For QPCs the scattering in¬ 
duced by disordeiP and the inhomogeneity of the edges— 
destroy the conductance quantization which can still be 
observecPS for mesoscopically smooth boundaried^of the 
constriction. As compared to QPCs, we find that an 
asymmetric narrowing, i.e. a contact between nanorib¬ 
bons of varied width, exhibits a very clear conductance 
quantization as a function of the Fermi energy, also when 
the edges are not smooth at the mesoscopic scale. We 
study quantum point contacts and find that for short 
constrictions there is a flat minimum of conductance for 
the tip above the constriction in agreement with the re¬ 
sults of RefP^l. Generally, for any tip potential we find 
that outside the graphene constrictions the variation of 
the conductance by the tip p oten tial is not pronounced 
as in the experiments of Refs.LUM3 p or longer constric¬ 
tions however a regular formation of the standing waves 
occurs which - according to the present study - should 
be resolved by the scanning gate microscopy. 


The resonant states localized at armchair-zigzag con¬ 
nections at ribbon constrictions formed by atomic 
step d 22 * 4 ^ induce a strong backscattering. These resonant 
states quench the conductance to zero at low energy. We 
find that the charge transport can be unblocked by the 
potential of the probe which neutralizes the scattering 
centers, with conductance maps forming halos similar 
to those found in the Coulomb blockade experiment of 
RefP^. The response for homogenous nanoribbons to the 
probe potential at low Fermi energies is determined by 
the existence of perfectly conducting channels (metallic 
armchair and zigzag ribbons), as well as by formation of 
local n-p junctions at the edges (zigzag ribbons). 


This paper is organized as follows: in Section II we 
outline the calculation method which is described in de¬ 
tail in the Appendix, Section III, IV and V describe the 
response of the nanoribbons, their narrowing and con¬ 
strictions, respectively (see Fig. [I]). Summary and con¬ 
clusions are given in Section VI. 



FIG. 1: Schematic drawing of the systems studied: nanorib¬ 
bon (a), asymmetric narrowing (b) and QPC (c). 


II. METHOD 


u=-l u =0 



FIG. 2: Schematic drawing of the nanoribbons in the input 
and output channels. In the Bloch waves that serve as the 
eigenstates of the channels u numbers the elementary cells, 
v the ions within the cell acc = 0.142 nm is the nearest 
neighbor distance, a = 0.246 nm is the lattice constant, and 
Ax is the period of the considered nanoribbon. 


We use the tight binding Hamiltonian for n electrons 
with the nearest neighbor hopping 

H = + c h ) + y ( r *)4 c i, (i) 

i*j} » 


and the hopping energy of t t j = —2.7 eV. For the de¬ 
scription of the transport at the Fermi level we pre¬ 
pared a tight-binding version of the quantum transmit¬ 
ting boundary method 37 (for the finite difference ap¬ 
proach see Ref.®). The studied system consists of an 
input and output channels and a scattering region. We 
assume that the channels are infinite and uniform, i.e. 
periodic in the atomic scale. The electron eigenstates in 
the nanoribbons are determined in the Bloch form, 

k u, v = X k v e ikuAx , (2) 


where k is the wave vector, u numbers the elementary 
cells of the ribbon (see Fig. [2]), x; — the atoms within the 
elementary cell, and Xv is the periodic function, assuming 
the same values within each elementary cell of the ribbon. 
The detailed procedure for determining the dispersion 
relation and Bloch functions for the channels is given in 
the Appendix. 

Generally, in the Hamiltonian eigenstates within the 
channels the scattering wave function is a linear combi¬ 
nation of the Bloch functions for a number of subbands 
at the Fermi level 


^u,v = + dlx k v -e ik - uA *). (3) 

i 


For the electron incident to the scattering region from the 
left in the input (left) lead far away from the scattering 
region the electron wave function takes the form of a 
superposition of incoming (k + ) and backscattered (&_) 
Bloch waves, 


ipo,v 


E ( I k\ I k l _. 

\CinX v ):> 

I 


( 4 ) 
















3 


where c\ n are the amplitudes of the incoming states and 
d\ n - of the outgoing ones. The sum in Eq. Q runs over 
the current carrying states since at a large distance from 
the scatterers the contribution of the evanescent modes 
can be neglected. In the right lead the wave function is 
a superposition of Bloch waves traveling from the scat¬ 
tering region, 


i’N,v = '£ C outXv + e< NAx . 


( 5 ) 


For solution of the scattering problem we assume that 
the electron is incident in a single subband l, c\ n = Sji. 
Then, the amplitudes d\ n and c l out are to be determined. 
At the boundaries between the leads and the scattering 
region which we cover within the computational box we 
require continuity of both of the wave function and its 
derivative. The numerical procedure for determination 
of the scattering amplitudes is given in the Appendix. 

Once the scattering amplitudes are determined we pro¬ 
ceed to calculation of transmission T and reflection R 
probabilities for th e zth incident mode according to the 
Landauer approac H 42 ! 43 !, T) = J2 tij and Ri = X) r ij > 

i j 


where tij = 


(4>k 


and rij = 


dL 


{4>k Y 


Wk+y 


are 


the probabilities of the electron transfer from incident 
subband i to j subband of the output and input chan¬ 
nels, respectively, and (fik+Y [(</>&_)*] is the probability 
current flux of the ith mode wave function traveling in 
the right [left] direction. 

The currents in the tight binding approach flow along 
the 7 r interatomic bonds, and the formula for the current 
flowing from atom m to atom n as derived^ from the 
Schrodinger equation is 


2 

jmn = , \tmnYmY' r> ) (d) 

where ipm is the wave function at the ?n-th ion. The flux 
is evaluated as 


<l> = '52'52jin„ ( 7 ) 

l ni 

where the sum runs over the atoms across the end of 
the channel l and their neighbors rq within the interior 
of the computational box. Finally, the conductance is 
evaluated from the Landauer formula as G = T). 

Since the formula involves sum over incident subbands, 
for the calculations it is sufficient to consider the solution 
of the Schrodinger equation for a single input channel 
only. Note, that our choice of the quantum transpar¬ 
ent boundary conditions implies working with the wave 
functions - whic h is an alternative to the Greens func¬ 
tion methocP2l22l, that deals with Hamiltonian operators 
mainly. The choice between the two is a matter of taste. 

The effective potential of the charged tip as seen by 
the Fermi level electrons is a result of the screening of its 


Coulomb potential by the electron gas. The effective po¬ 
tential as obtained by the Schodinger-Poisson modeling^ 
can be approximated by a Lorentz function 

V{X ’ V)= 1 + ((* - x t y + (y - Vt) 2 ) /cP ’ (8) 

where x t ,yt are the coordinates of the tip position, d 
- the width of the effective tip potential and V t - its 
height. The height of the tip potential is determined by 
the charge accumulated by the tip, and d is of the range 
of the tip - electron gas distance™ A discussion of the 
results for varied Vt and d parameters is provided below. 


III. CONDUCTANCE MAPPING FOR 
NANORIBBONS 

We consider first the zigzag ribbons - labeled by ZZ in 
the Figures, and next proceed to armchair: semiconduct¬ 
ing (AS) and metallic (AM) ribbons. 


A. Zigzag ribbons 

Figure [3](a) shows the dispersion relation for a zigzag 
nanoribbon with 102 atoms across the channel. The 
zigzag edge does not couple the states of the K and 
K' valleys of graphene which are thus present in the 
dispersion relation (K,K' = ±|^). The dispersion 
relation is given in Figs. |3][a-c). For Fermi energy 
Ep = 0.44 eV there are five eigenstates with current 
flowing to the right (kp , ky , k 3 > , k 2+ , ^ 3 +) and to the 
left [ki _, k 2 '_, k 3 '_, k 2 _, k 3 _). 

Let us consider the ribbons conductance response to 
the tip potential for a low energy of the incident elec¬ 
tron. Figures |4][a-b) show the transfer probability as a 
function of the tip position across the ribbon for two 
values of Ep corresponding to the lowest subband trans¬ 
port (electron incoming with wave vector kp , that can 
be backscattered only to k\_ - see Fig. [3](b,c)). For both 
the considered Fermi energies the conductance gets low 
when the tip is located near the edge, in spite of the fact 
that the current for the unperturbed zigzag ribbon van¬ 
ishes at the edges [FigureQd-e)]: at the zigzag edges only 
the ions of a single sublattice are occupied in the Hamil¬ 
tonian eigenstates [Figure |4][g-h)] and the current flows 
between the ions belonging to different sublatticeJ^Sl. For 
the central position of the tip we have either T = 1 for 
larger Ep = 0.13 eV [Fig. |4][b)] or a rapid variation of T 
between 0 and 1 for the lower Ep = 0.07 eV [Fig. |4][a)]. 
In the plot of G as a function of the tip position and the 
Fermi energy given in Fig. [5] we find a wedge-shaped flat 
region with T = 1. This region is related to the perfectly 
conducting channel 17 18 in the zigzag ribbon. Perturba¬ 
tion by the tip - even if large - does not induce scattering 
between the K and K' valleys as long as the tip potential 
is slowly varying at the atomic scale. Since the tip po¬ 
tential has a long range character there is no intervalley 
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scattering, and in consequence in the lowest subband 
transport conditions - the backscattering is absent, un¬ 
less the tip is located near the edge of the ribbon [cf. Fig. 
[6]. The tip placed in the center of the ribbon - is not 
exactly transparent for the electron flow, since the cur¬ 
rent makes its way between the edge [Figure 0 a)], but 
anyway no conductance response is observed. For the tip 
near the edge, the current circulates around the defect 
before it is backscattered (see Fig. 0b) - the illustrated 
case corresponds to T = 0). The edge mediates in the 
intervalley scattering provided that its potential is raised 
to the Fermi energy® by e.g. the tip inducing formation 
of a local n-p junction. In Fig. [5] we marked with the 
black line the position of the tip for which its potential at 
the lower ribbon edge is raised to the Fermi energy, with 
a perfect agreement with the boundaries of T = 1 wedge 
region. The transfer probability as a function of the tip 
position plotted in Fig. 0a) for low Ep exhibits a rapid 
variation between 0 and 1 which corresponds to the case 
when the tip potential at both edges exceeds Ep, hence 
multiple intervalley scattering events occur, resulting in 
either T = 0 or T = 1. 


The T = 1 wedge region continues smoothly into the 
regime of a few transport subbands at the Fermi level 
[see Fig. [5], For larger energies, in the right (K) valley 
of the dispersion relation [Fig. 0a, c)] the number of 
modes carrying the current to the right is by 1 larger 
than the number of modes carrying the current to the 
left [Fig. 0c)]. Hence, even for a strong but long-range 
pert urbati on the ribbon carries current at least in a single 
modc 1 ®^ unless the tip raises the edge potential to the 
Fermi level - as it is still the case in Fig. [5ja,b) also 
beyond the single subbands transport conditions. 


In Figure[8]we plotted the backscattering probabilities 
for the parameters of the tip of Fig. 0a) (large and 
wide tip potential V = 0.5 eV, d = 40 nm) for the Ep 
range with three subbands or six wave vectors at the 
Fermi energy [Fig. 0b- c)]. The intervalley scattering 
in the lowest subband (of probability Rp\ corresponding 
to the change of the wave vector from k' 1+ to - see 
Fig. 0b, c)) occurs only for the tip near the edge of the 
ribbon [Fig. 0a)], but when it does, it is nearly complete 
(Rpi — 1). The other intervalley scattering paths Ry 2 
[Fig. 0 c)j, i? 2'2 [Fig. 0e)], f? 22 ' [Fig. 0f)] occur also 
for the tip near the edge but with smaller probabilities. 
The effective backscattering channels for the tip in the 
center of the ribbon are R 2 ’ 2 : = R 22 [Fig. [8][d)J that 
appear within the same parabolic subbands [Fig. f3Tb,c)] . 
The backscattering by the tip in the k 2+ —> k' 2 _ channel 
governs the current and scattering density pattern that 
is displayed in Fig. 0c-f). The oscillations at the left to 
the tip are due to the interference between k 2 ' + and k 2 ' 


eigenstates and correspond to the period of | fc t _ kl \ , 

which is equal to 103A and 58A for Ep = 0.3 and 0.4 eV, 
respectively. The conductance of the ribbon raises from 
1 to 3 for Ep above 0.2 eV [see the dashed lines in Fig. 
[5j with backscattering by the tip in the 2' + —► 2'_ channel 


disappearing at higher Fermi energies, [cf. Fig. UN- 


B. Armchair ribbons 

The armchair edges of the ribbons strongly mix the val¬ 
leys and the zero point energy appears only for a single 
wave vector in the dispersion relation for metallic chan¬ 
nels [Fig. 0d)]or a single extremum per energy band is 
observed [Fig. [31(g)] for the semiconducting channel. 

For the semiconducting armchair ribbon we find that 
the tip strictly blocks the current for any low Ep [region 
T’ in Fig. 0b)] while the metallic ribbon ignores the 
presence of the tip in the single-subband transport con¬ 
ditions [region T in Fig. 0a)] . In contrast to the case of 
the zigzag ribbon, the conductance is a smooth function 
of both Ep and the tip position. The low-energy results 
are given in Fig. [4] with T = 1 for the metallic ribbon 
for both Ep = 0.07 eV and Ep — 0.13 eV, irrelevant of 
the tip position. The metallic armchair ribbon possesses 
a perfectly conducting channel due to the pseudovalley 
symmetry®, conserved at low energy, hence the absence 
of backscattering in the lowest subband transport. For 
the semiconducting armchair ribbon the tip blocks the 
current for any tip position above the ribbon [Ep = 0.07 
eV in Fig. 0a)] or reduces it strongly, particularly near 
the center of the ribbon [Ep =0.1 eV in Fig. 0b)]. In 
general, the most effective backscattering channel is the 
one within the same subband [k 2+ —> k2~, k 3+ ->■ fc 3 _] 
and for the tip above the center of the ribbon. 

At higher energies, for the semiconducting ribbon, 
only the backscattering within the same subband is non- 
negligible [Fig. |l0|(d-f)], including the non-zero chan¬ 
nel in the lowest subband. For the metallic ribbon the 
backscattering within the lowest subband is absent due 
to the pseudovalley symmetry®. Scattering in the same 
subband [Fig. [Toj(b-c)] is the strongest as in the semicon¬ 
ducting case and it is most effective for the tip above the 
center of the ribbon [Fig. [Tojb-f )]. For the metallic rib¬ 
bon a non-zero scattering between the subbands is also 
found [Fig. [To|[a)] but only for the tip near the edges of 
the ribbons. In Fig. 0c) we plotted a cross section of 
Fig. 0c ) and Fig. 0] for Fermi energies in multiple sub¬ 
band transport regime. For higher Fermi energies, the 
tip can no longer form the n-p junction at the edge, and 
the characteristic backscattering of the zigzag ribbon is 
not observed. 

Summarizing, for higher Fermi energies the overall re¬ 
sponse of the ribbon conductance to the tip position is 
qualitatively similar for all ribbon types, i.e. it is smooth 
and maximal at the center of the ribbon. The T(y t ) de¬ 
pendence can be approximately put in a (yt — D/ 2) Q form 
with a that strongly varies with the energy for any type 
of the edge, D being the width of the ribbon. For higher 
Fermi energies the response of the ribbon to the tip scan¬ 
ning its surface is qualitatively independent of the type 
of the boundary, in spite of the fact that the current and 
density distributions are distinctly different [Fig. 0f,i)]- 
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FIG. 3: (a) Dispersion relation for a zigzag nanoribbon 102 atoms wide. (b,c) Enlarged fragments marked by rectangles in 
(a). The arrows indicate the current direction, (d-f) same as (a-c) only for a metallic armchair ribbon with 92 atoms across 
the channel, (g-h) same as (a-c) only for a semiconducting armchair ribbon with 93 atoms across the channel. Arrows in 
(b,c,e,f,g,h) indicate the direction of the current flow. Here and in the other figures ZZ, AM, and AS stand for zigzag, armchair 
metallic and armchair semiconducting ribbons. 



FIG. 4: (a-c) The transfer probability as a function of the tip position for Vt = 0.2 eV and d = 40 a. (d-f) The density current 

averaged over neighbor cells, (g-i) The probability density for ribbons without the tip. In the left column (a,d,g) Ef = 0.7 eV 
and in the middle (b,e) Ef = 0.13 eV for zigzag and armchair metallic and (h) Ef = 0.1 eV for armchair semiconducting. Red 
lines correspond to a zigzag ribbon with 102 atoms across the channel, green (blue) to the metallic (semiconducting) armchair 
ribbons with 92 (93) atoms across the channel. In (c,f,i) results for higher energies are given. For the corresponding dispersion 
relation see Fig. |3|a-c) for the zigzag ribbon, and Fig. §d -f) for metallic and Fig. §g- i) for semiconducting nanoribbons. 


IV. NARROWING OF THE RIBBON 


As an intermediate step towards description of QPCs 
let us consider the channel of varied width that is de¬ 
picted in the inset to Fig. 0 The channels at both sides 
of connection (A,B and E,F in Fig. 11) as well as the 
linking edges ( C, D in Fig. 11) are assumed of the same 
type (zigzag or armchair). Short segments of varied type 
of the edge appear at the connections between the hor¬ 


izontal and slanted ends of the ribbon (see the inset in 

Fig. pTi)). 


Figure 0 shows the conductance of the contact as a 
function of the Fermi energy for both zigzag and armchair 
edges. The conductance exhibits a very neat quantization 
as a function of the Fermi energy. A pronounced quan¬ 
tization of conductances j s generally difficult to obtain 
in graphene quantum point contacts (see Introduction). 


According to the results of Fig. 11 a channel contain- 
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FIG. 5: Conductance of the zigzag ribbon with dispersion relation of Fig. [3] as a function of the position of the tip and the 
Fermi energy. Parameters of the tip Vt = 0.5 eV, d = 40a (a), Vt = 0.5 eV, d = 20a (b), Vt = 0.2 eV, d = 40a (c), Vt = 0.1 eV, 
d = 40a (d-f). The dashed lines separate regions of a varied number of subbands at the Fermi level - given above the upper 
axis. The thick black line in (a-d) corresponds to the tip positions for which the potential at the edge is raised to the Fermi 
energy. 
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FIG. 6: Scattering probability density for Vt = 0.2 eV and 
d = 40a at Ef = 130 meV and the tip near the lower edge 
of the ribbon (cross). The circle corresponds to potential 
introduced by the tip equal to Ef- 
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FIG. 7: (a,b) Amplitude of the current for the zigzag ribbon 
for V t = 0.2 eV, d = 40 a for energy Ef = 130 meV in the 
lowest subband transport and the tip in the center of the 
ribbon (a) - for which T = 1 and near the lower edge (b) 
of the ribbon [T = 0]. (c-f) results for three subbands at 
the Fermi level for the amplitude of the current (c,e) and the 
scattering density (d,f) for Ef = 0.3 eV and Ef = 0.4 eV. 
The summed transfer probability is 1 at (c,d) and 3 at (e,f). 



FIG. 9: Conductance versus the Ef and tip position for a 
metallic (a) and semiconducting (b) nanoribbons for Vt = 0.2 
eV, and d = 40 a. The dashed lines separate the regions with 
varied number of subbands at the Fermi level. 


ing a narrowing is a good candidate for observation of 
conductance quantization in graphene systems. 


A. Zigzag system 


FIG. 8: Intersubband backscattering probabilities for pa¬ 
rameters of Fig. § a) a) Ri'i, b) c) Ri> 2 , d) R 22 , e) 

I? 2 ' 2 , f) R 22 ' for Vt = 0.5 eV, d = 40a. 


Let us consider the narrowing of the ribbon with zigzag 
edges. At low Fermi energy of 50 meV the connection of 
zigzag edges has a very low conductance of T = 0.006 [see 
the orange line in Fig. [IT]. The ribbons at both sides of 
connection contain a perfectly conducting channel, but 
the inhomogeneity of the width leads to a very strong 
backscattering. At the connection of the zigzag edges a 
short sequence of an armchair form appears [see the arrow 
in Fig. i)], which opens the intervalley backscattering 
channel. The tip potential above any of the two points 
closes the intervalley backscattering channel, allowing for 
the electron passage to the wider ribbon [Fig. |l2|(a,b,i,j)]. 
When the tip is close to the zigzag edges, it produces a 
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FIG. 10: Intersubband scattering probabilities for armchair 
metallic (a,b,c) and semiconducting (d,e,f) nanoribbons for 
Vt = 0.2 eV, d = 40a. The number of subbands at the Fermi 
level is given above the plots. For the metallic ribbon (a-c) 
Tin is zero. 



FIG. 11: Conductance at the contact between channels of 
varied widths (inset). The channels at the left (A,B) and right 
(E,F) hand side of connection and the edge at the variation 
of the width (C,D) are of the same type : zigzag or armchair. 
The lines show the conductance for zigzag, armchair metallic 
and semiconducting (offset by 1 in T ) armchair edges of the 
channels. The number of atoms across the narrow (wide) 
channels is 188 (564) for the zigzag system, 161 (363) for 
metallic armchair channels and 163 (365) for semiconducting 
armchair edges. The dots on the lines indicate the points for 
which the SGM scans of Figs. 1 12| and [l4| were calculated. 


rapid variation of conductance which is related to forma¬ 
tion of the n-p connection at the edge as discussed above. 
Note, that this applies not only for the narrower ribbon, 
but also at the widening and at the wider ribbon. Thus 
the probe resolves the zigzag edges in the entire system. 
The strong conductance response near the scattering cen¬ 
ters localized at the edges are the counterparts of the 
conductance halos found in SGM mapping of graphene 
constrictionsFor the tip above the center of the nar¬ 
rower ribbon in Fig. [l2|)b) the potential at the edge stays 
below the Fermi level and the presence of the tip is not 



FIG. 12: Conductance maps for connection of zigzag channels 
of varied widths for Vt = 0.1 eV and d = 20a (a-d) and Vt = 
0.2 eV and d = 40a (e-h). Panels (i-1) show the scattering 
probability density for the electron incident from the left in 
the absence of the tip. By arrows in (k.l) we indicated the 
period of the probability density oscillation which corresponds 
closely to the values expected for the backscattering within 
the last subband (see text). The inset to Fig. 12(i) shows the 
atomic structure of the turn of the edge with a short armchair 
segment that induces backscattering at low Fermi energy. 


resolved in the map which is flat in the interior of the 
narrow ribbon [Fig. [T2](c,f )]. The insensitiveness of con¬ 
ductance to the tip position in the center of the narrow 
channel and the rapid variation at the zigzag edge - are 
consistent with the properties of the zigzag ribbons for 
low Fermi energies discussed above. 

For higher energies the effects of the edges are no longer 
present in the conductance maps [Fig. 12 c,d,g,h)]. In¬ 
stead we observe - particularly for the lower potential at 
the tip [Fig. [l2|[c,d)] an appearance of a periodic oscilla¬ 
tion of conductance. For instance the T maxima in the 
SGM maps are spaced by 202A for Ep = 0.160 eV and 
coincide with the minima of i? 2 ' 2 ' = R 22 backscatter¬ 
ing. The T maxima are spaced by a characteristic wave¬ 
length for the intersubband scattering Ax = ^ _ k ^ | . 


The wave vectors are equal to /c 2 + = 0.6825^ and 
& 2 _ = —0.6582^ for Ep = 0.160 eV, producing spac¬ 
ing Ax equal to 202.3A, in a perfect agreement with the 
oscillation period found in the SGM map. For larger Vt 
and wider d, a deviation from this perfect agreement is 
found due to variation in the wavelengths induced by the 
external potential of the tip and backscattering involving 
a larger number of subbands. In particular in Fig. |12[ g) 
for V t = 0.2 eV the transport is reduced to the lowest 
subband (T = 1) for the tip above the channel and the 
oscillations disappear in the SGM map which takes the 
form similar to the one found for lower Fermi energy - 
Fig. [l2](a,b). 

The oscillations of the scattering density observed 
above for the uniform ribbons [see Fig. [ 7 ] resulted from 
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the presence of the tip but were not resolved by the con¬ 
ductance mapping. For the non-uniform ribbons the os¬ 
cillations appear in the scattering wave functions due to 
the inhomogeneity of the channels already in the absence 
of the probe [Fig. [l2|(k,l)] and according to the present 
results they can be resolved by the scanning gate mi¬ 
croscopy. 

The number of subbands that appear at the Fermi 
level increases with the energy and several backscatter- 
ing channels appear. For Ep = 280 meV [Fig. [l2^d)] 
the spacing between the extrema of T map are ~ 108A 


The backscattering probabilities are displayed in Fig. 13 


The SGM map T for E = 280 meV inside the channel 
[see Fig. [l2|[d)] is very well correlated to spatial variation 
of backscattering R 3 ' 3 > = R 33 . The oscillation period re¬ 
sulting from the difference of the wave vectors is 107.23 A 
with a perfect agreement with the spacing of T extrema 
in the SGM map. The intervalley scattering [i? 3 / 3 in Fig. 
[M b )] - gives only an enhanced backscattering at the 
exit from the thinner channel with no pronounced oscil¬ 
lation within. We find as a general rule, that - at least 
for a moderate tip potential [V t = 0.1 eV, d = 20a] 
the oscillations of conductance are observed as a func¬ 
tion of the tip position inside the zigzag thinner channel 
and have periodicity corresponding to the backscattering 
within the highest subband. 


B. Armchair systems 


Figure [14] displays the T maps for the armchair edges of 
the narrowing [see the inset to Fig. |TT] . At low energy the 
conductance of both metallic and semiconducting ribbons 
is nearly zero [Fig. 11 . Figure |l4|[a,d) indicates that 
similarly as for the zigzag system at low energy [Fig. 
|l2|[a,b)] the conductance is raised when the tip is located 
at the turns of the edges: at connection of the slanted 
edge to the thinner (for the metallic ribbon Fig. [l4](a)) or 
wider (for semiconducting ribbons Fig. |14[ d)) armchair 
channel. 

For Ep = 90 meV and the thin semiconducting ribbon 
we have already two subbands at the Fermi level [Fig. 11 
and the T map contains a clear oscillation [Fig. 
of periodicity characteristic to the backscattering in the 
uppermost subband. For the metallic [Fig. rib¬ 

bon at this energy the thinner ribbon is still in the single 
subband transport conditions, no oscillations occur - ex¬ 
actly as for the zigzag ribbon for low Fermi energy [see 
Fig. 12 ]a-b,e-f)[. For armchair systems at Ep = 90 meV 
no edge effects are observed and the backscattering oc¬ 
curs only for the tip near the exit from the thin channel. 
The oscillations for the metallic ribbon occur when the 


next subbands enter below the Fermi energy [Fig. 11 


see Fig. |l4[c), with the periodicity characteristic gener¬ 
ally to backscattering within the highest subband with 
the period of | fc | for l being the index of the last 
subband. There are exceptions to this rule, in particular 
when the subband is very close to the Fermi energy the 


backscattering of the subband can be nearly complete, 
then the one which is observed in the oscillation is not 
the highest but the next lower one. 


0.13R 33 - R 3'3'0 0.27 
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FIG. 13: Intrasubband (a) and intervalley (b) backscattering 
for the zigzag narrowing with the Fermi energy of 280 meV 
for V t = 0.1 eV, d = 20a. 
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FIG. 14: T maps for a narrowing with armchair edges: metal¬ 
lic (a-c) and semiconducting (d-f) for Vt = 0.1 eV and 
d = 20 a. The insets in a) and d) display the fragments of 
the edges near the exit from the narrower channel and the 
entrance to the wider channel, respectively where a short seg¬ 
ment of zigzag edge appears, leading to the strong backscat¬ 
tering which is neutralized by the probe potential as in Fig. 
QJ above. 


V. QUANTUM POINT CONTACTS 
A. Long narrowing of QPC 

We are now ready to proceed to description of SGM 
maps for quantum point contacts formed within the 
graphene ribbons. We first consider a system based on 
two narrowings of the channels as described above, in the 
form given in the inset to Fig. [l5]with a finite length L 
of the narrower ribbon. We assume that the narrower 
ribbon has a width d„ — 20 nm and consider L between 
32 nm and 58 nm, for the transverse size of the wider 
ribbon between 44 nm and 60 nm. The conductance of 
all the considered ribbons - zigzag and armchair - both 
semiconducting and metallic - drops to zero for low Ep. 
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FIG. 15: Transfer probability as a function of Fermi energy 
for QPC with long constriction (inset). The dashed lines rep¬ 
resent the transfer probability in graphene ribbon of the same 
width as the constriction. The conductance for armchair rib¬ 
bons is shifted by 1 for the metallic system and 2 for semi¬ 
conducting system. The lengths indicated in the inset are 
L = 577 A, d n = 199.26 A and d w = 447.72 A for the arm¬ 
chair semiconducting ribbon, L = 551 A, d n = 196.8 A and 
d w = 445.26 A for the armchair metallic ribbon and L = 320 
A, d n = 198.8 A and d w = 599.24 A for the zigzag ribbon. 


The inhomogeneity of the ribbon leads in this way to ap¬ 
pearance of the transport gapPS for any type of the edge 
of the wider ribbon. 

The conductance of the QPC with zigzag edges is dis¬ 
played in Fig. [15] with the solid orange line, while the 
dashed orange line indicates the conductance of the nar¬ 
row part of the system (the number of subbands carrying 
the current to the right at the Fermi level). The conduc¬ 
tance of the QPC is generally lower than the conductance 
of the narrower part with exceptions of a series of con¬ 
ductance peaks. In Fig. M we plotted the scattering 
probability density [Fig. 16 ai-ap)] and the SGM maps 
[Fig. 16 [aa-ah)[ for the zigzag QPC and a number of 


chosen Fermi energies marked in Fig. [15] with points. In 
the T maps we find a pronounced variation of conduc¬ 
tance maps within the narrower part of the ribbon. Out¬ 
side the constriction the conductance maps are almost 
flat. No conductance variation of visible in the scale of 
the variations within the narrower channel. For the low- 
energy peak Ep = 76 meV - within a single subband of 
the narrow channel - the conductance map of the system 
drops when the probe is located within the constriction 
[Fig. [I6|(aa)]. In the scattering densities and conduc¬ 
tance maps we find signatures of backscattering at the 
turns of the edges as discussed above for the asymmetric 
narrowing of the channel [Fig. |16[ aa)]. For higher ener¬ 
gies [Fig. 16 ab-ah)] the corners correspond to high val¬ 
ues of the probability densities but with no counterparts 
in the SGM map - again as in the asymmetric narrowing 
of Fig. [12} 

For higher Fermi energies, within each segment of Ep 
corresponding to a fixed number of subbands of the nar¬ 
row channel (dashed orange line in Fig. 15 we observe 
that 1) at the first maximum of T(Ep) at the low en¬ 
ergy side the conductance map contains a single mini¬ 


mum along the channel, 2) each subsequent T peak that 
is higher at the Ep scale corresponds to conductance 
map with a number of minima increased by 1. A rep¬ 
resentative example is given in Fig. [l6|[ab-ad) for the 
T maps for the first, second and the third T maxima of 
conductance [see the numbers near the orange line in Fig. 


15 . The conductance resonances correspond to an inte¬ 


ger number of the scattering wavelengths for the highest 
subband within the constriction with maxima observed 
in the scattering probability density [cf. Fig. 16[aj-al)]. 
Once another subband falls below the Fermi energy (for 
instance above Ep = 0.2 eV) another sequence of conduc¬ 
tance peaks is formed with the number of maxima within 
the constriction starting from 1 again (see the peaks la¬ 
beled by primed numbers and maps of Fig. 16[ae,af)]. 
Beyond the resonances this regular correspondence be¬ 
tween the scattering density and the conductance map is 
lost [cf. Fig. 16 ag,ah) for energies Ep = 146 meV and 
Ep = 177 meVA 

The properties of the QPCs with armchair edges and 
the width corresponding to a semiconducting dispersion 
relation turn out to be similar to those of the zigzag sys¬ 
tem, with pronounced T(Ep) peaks [cf. blue solid line in 
Fig. 15 corresponding to a determined number of max¬ 
ima of the probability density correlated with the con¬ 
ductance maps. For the metallic armchair system T(Ep ) 
dependence is more complex [red line in 15 . In particular 
we find that the peaks corresponding to a fixed number 
of extrema of the probability density and SGM maps ap¬ 
pear in pairs [cf. Fig. 16[l,m) for Ep = 93.5 meV and 
Ep = 96.3 meV for a single extremum, and Fig. 16 [n,o) 


for Ep = 104.9 meV and Ep = 109.3 meV for the double 
extremum, etc.]. 


The difference in the properties of semiconducting and 
metallic armchair systems has its source in the details 
of the dispersion relation. When a subsequent subband 
for the electron transport is opened within the narrow 
channel, the subbands of the metallic ribbon are nearly 
degenerate [see Fig. [3](d) for low k\ while for the arm¬ 
chair ribbon the subbands are split near zero wave vector 
[see Fig. 3j[g)] . Thus, for the metallic ribbon we have 
two close highest subbands of similar Fermi wavelengths, 
inducing formation of separate but close resonances for 
both. The neighborhood of the wavelengths enhances 
the intersubband scattering leading to a complex form of 
T(Ep) at higher energies. 

Concluding, for QPCs with a long constriction the con¬ 
ductance peaks of T(Ep) dependence are related to reso¬ 
nances of the highest subband, and the resonant density 
within the constriction is resolved by the SGM conduc¬ 
tance maps. The tip potential enhances the backscatter¬ 
ing within the highest - nearly parabolic - subband. The 
tip-induced backscattering is most effective in the region 
where the scattering probability density is maximal. The 
zigzag and semiconducting narrowing display a distinct 
regularity in the subsequent number of peaks and the 
form of the SGM map. For the metallic armchair ribbon 
the T[Ep) dependence is more complex due to the near 
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FIG. 16: T maps for a long QPC with zigzag edges (a-d) and scattering densities (e-h) at the Fermi energies corresponding to 
the maxima of transmission probability. Parameters of the tip: V t = 0.1 eV, d = 20 a. 


degeneracy of the highest subbands near their bottom 
and pairs of resonances with similar conductance maps 
are found. For each system, outside the constriction the 
conductance map is almost flat. 


B. Short constriction 



FIG. 17: Transfer probability as a function of Fermi energy 
for wedge-shaped QPC (inset). The widths indicated in the 
inset are d n = 58 A, d w = 598.8 A. 


Let us consider the QPC with a short constriction. We 
studied first a wedge shaped constriction of a zigzag rib¬ 
bon depicted in the inset to Fig. [17] with the ribbon 
narrowing from about 60 nm to about 6 nm. The con¬ 


ductance of the QPC is displayed in the main panel of 
Fig. |17| The short constriction cannot support any lo¬ 
calized resonances in the form discussed above, and the 
corresponding series of T(Ep ) peaks are not observed. 
For any Fermi energy a minimum is found in the SGM 
map within the constriction and only a low variation of 
conductance outside the constriction is found - see the 
maps displayed in Fig. [l8j Some variation outside the 
constriction is only observed for low Ep [Fig. [l8[a)] and 
it disappears at higher energies in the scale of the varia¬ 
tion in the narrow part. [Fig. 18 b,c)]. 

The results found for this short, narrow and abrupt 
QPC are reproduced also for a large smooth at the 
atomic scale constriction of a cosine profile^. The system 
exhibit^ a detectable smooth steps of conductanceP-0 
also when the constriction remains quite wide. The 
T(Ep) result is presented in Fig. 19 where we considered 
wide (d w = 100 nm to 150 nm) zigzag and semiconduct¬ 
ing armchair edges that at the narrowest part have width 
of d n = 33 nm to 50 nm. 

The conductance maps for the zigzag ribbon are dis¬ 
played in Fig. [20j In the scattering probability den¬ 
sity maps we find maxima along the constriction [Fig. 
[20[e-h)]. The atomic steps forming constrictions contain 
multiple of zigzag-armchair short segments of the edges 
[Fig. [20[i)j - which as discussed above - lead to an effec¬ 
tive backscattering. At low energy [Ep = 25 meV in Fig. 
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FIG. 18: T maps for a wedge-shaped QPC with zigzag edges (a-c) and scattering densities (d-f). Parameters of the tip: 14 = 0.1 
eV, d = 20 a. 


[20^a)] the tip neutralizes the backscattering when located 
above these parts of the edges - in consistence with the 
results for asymmetric narrowing [cf. Fig. [l2{a,b)]. At 
higher energies [cf. Fig. [20|[b-d)] we are left with a mini¬ 
mum of conductance at the center of the constriction and 
a low variation of conductance in the outside. 



FIG. 19: Transfer probability as a function of Fermi energy 
for large QPC (inset). The widths indicated in the inset are 
d n = 497 A, d w = 1498.1 A for zigzag ribbons and d n = 332.1 
A, d w = 998.76 A for armchair ribbons. 


VI. SUMMARY AND CONCLUSIONS 

We have performed an analysis of the conductance re¬ 
sponse of graphene nanoribbons and their narrowing to a 
perturbation introduced by potential of a scanning probe. 
The study was based on the tight binding implementa¬ 
tion of the quantum transmitting boundary method. 

We find that at low Fermi energies the response of the 
zigzag ribbons to the scanning probe occurs only for the 
tip above the edges of the channels due to formation of 
the local n-p junctions allowing for a very effective inter¬ 
valley backscattering and producing rapidly varying con¬ 


ductance maps. For the probe above the interior of the 
zigzag ribbon the conductance does not react to the per¬ 
turbation. The finding is due to existence of the perfectly 
conducting channel and is also observed for metallic arm¬ 
chair ribbon irrespective of the position of the probe. 
For higher Fermi energies the response of the zizgag and 
armchair nanoribbons is similar producing smooth con¬ 
ductance maps as functions of the tip position with the 
strongest response in the center of the ribbon. 

For the asymmetric narrowing of the ribbon - which 
exhibits a very distinct conductance quantization - the 
conductance maps at low energy resolve the zigzag edges 
and scattering centers near the zigzag-armchair segments 
which necessarily appear in atomic steps that form the 
narrowing. For higher Fermi energy an oscillation of 
the scattering probability density is observed, which in 
general corresponds to backscattering within the high¬ 
est subband of the narrower channel. According to the 
present results the oscillation can be resolved in the maps 
of conductance response to the probe. 

For quantum point contacts similar oscillations of con¬ 
ductance maps are observed for longer constrictions that 
are connected with resonant peaks of T(Ep) dependence. 
Formation of conductance peaks have a regular reentrant 
character in the energy with subsequent subbands ap¬ 
pearing at the Fermi level. The peaks correspond to 
scattering probability density and conductance map with 
a determined number of extrema along the constriction. 
Shorter constrictions that do not support the resonances 
within the narrowing produce a smooth T(Ep) depen¬ 
dence with SGM maps producing a minima of conduc¬ 
tance within the contrictions. For the asymmetric nar¬ 
rowing - at the wider part of the ribbon and for QPC 
outside the constriction no conductance response of the 
system is observed in the response maps. 

















12 


E f =25 meV 


E f = 60meV 


E F =85 meV E F =118 meV 





FIG. 20: T maps (a-d) for a large QPC with mesoscopically smooth edges within a zigzag nanoribbon and the corresponding 
scattering densities (e-h) in the absence of the probe. Parameters of the tip potential: V* = 0.1 eV, d = 20 a. (i) - enlarged 
fragment of the narrowing indicating armchair segments within the zigzag edge that are mainly responsible for the scattering 
at low Fermi energy and that are neutralized by the tip potential [see (a)]. 
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VII. APPENDIX 


A. Channel eigenstates 

The procedure for determining the dispersion rela¬ 
tion and the Hamiltonian eigenstates for the nanoribbon 
channels is described here. For an infinite, uniform lead 
the tight-binding Hamiltonian can be put in a matrix 47 


where matrix H„ concerns a single elementary cell u [see 
Fig. 1 and B u describes the connection between cells 
u and u + 1. The H„ and B u matrices have dimensions 
nxn where n is number of atomic orbitals within a single 
elementary cell. For ideally periodic nanoribbons forming 
the channels the matrices H u , B n are the same for all 
elementary cells, and we skip the u index in the following. 
The Hamiltonian eigenfunction inside the channel can be 
written in a form divided on separate cells 


( : \ 


■0 = 


4>u-i 

4>u 

V’u+l 


V ; 7 


( 10 ) 


where 0 U is a vector of size n. The wavefunction satisfies, 

- B0 U _! + (El - H)W + Bt0 u+1 = 0. (11) 

Based on the Bloch form of the wave function [Eq. [2] we 
make the following subsitution: 


fpu-1 = X, V’u = Ax, ip u+ 1 = A 2 x, 


for which 


-Bx + A(£I-H)x + A 2 B+x = 0. (12) 


With 


V = Ax 


the eigenproblem can be put in form 


( o I )_ x ( I 0 
[-B EI-H J A { 0 Bt 



(13) 


The generalized (2n x 2n) eigenproblem has 2 n solutions: 
n left-going and n right-going modes - propagating or 
decaying. 4 ' For evanescent modes it is straighforward 
to identify right and left going modes. For right-going 
evanescent modes the eigenvalue satisfies |A + ,„| < 1 and 
for left-going evanescent modes |A_, n | > 1. The propa¬ 
gating modes - for the Bloch waves of Eq. (pi) have the 
form A±,„ = exp(ikAx), with a real fc, hence~fA± jn = 1|. 
For a given E we look for the values of |A| = 1 and de¬ 
termine the corresponding wave vectors and the periodic 
functions X- Finally, we evaluate the current flux (p| and 
determine the direction of propagation the current flux 
is positive (negative) for right (left) going modes. 


B. Determination of the scattering amplitudes 



H = 


0 

0 


0 0 

H u -i Btr 0 

B u _i H u Bt 

0 B. u H u+1 

0 


\ 


0 

0 




The scattering amplitudes of Eq. Q and Eq. (p]) are 
determined in the following manner. We evaluate the 
coefficients d\ n multiplying (pj) by a complex conjugate 

(9) ofx 0( : 

<X fci - IV’o) = E( c »(x fc - lx fe +) + d\ n (x k - |x fe -». (14) 


V 


o 
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Here we denote the inner product in discrete form 
(A\B) = Defining vector A with A;/ = 

j 

<X fe -|V>o) and matrices = ( X k -\x k +), S;/,; = 


( X k -\x k ~) Eq. (14) can be written in a matrix form 


From the asymptotic conditions for the wave functions 
( 4|5 ) we derive the boundary conditions for the algebraic 
set of equations for c l out and d l iv . Let us consider deriva¬ 
tive of the wave function [Eq. Q] at the left end of the 
computational box 


A = Be, 


Sd, 


(15) 


which implies 

l' I'd 

d in = ~ C in2- ( 16 ) 

l’ 


1 

Aa; 


(ipo,v 


i>-i, v )=J2( c inxZ + 4z( 


Ax 

= E(°inXv + ^A kl+ + d l m Xv~ ^V). 


d‘ n Xv-^ x (l-e- ik - Ax )) = 


Using (16) we transform this equation to obtain ip^\y. 


V’-l.U 







(17) 


Similarly, for the other end of the computational box we multiply equation (|5j) by the complex conjugate of 
(X k +e ik + NAx ): 


( X k +e ik + NAx \ ijj N } = ^(x fc +|x' s +)e i(fe +- fe + )iVAx C L = E S J' 

l 


l^OUt 5 


(18) 


hence 


cl out = Y,&~ 1 )iAx <e<NAx \^n) 


The derivative at the end of the computational box is 

^z(i’N+ 1 ,„ - M =Y, cl out^ + e ik + NAx ^-(e i ^ Ax - 1 ) 


(19) 


( 20 ) 


Then, using (19) we obtain: 


1pN+ 1 ,d = 1pN,v + / , 

1,1' 


(s’-'Mx 


k + e ik + NAx \l/j N )Xv + 


e ik l + NAx ^ e i(k l + )Ax _ -g 


( 21 ) 


The boundary conditions given by and ( plj ) are used as expressions for the wave functions outside of the com¬ 
putational box in the Sc hrodinger equat ion Hip = Eip. Having calculated the wavefunction ip, we can calculate the 


coefficients d\ n from (16) and c l out from (fid)]) . The results are next used for evaluation of the scattering probabilities 


and conductance in the Landauer approach. 










